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Abstract 

A recent study of nonextensive phase transitions in nuclei and nuclear clusters needs a probability 
model compatible with the appropriate Hamiltonian. For magnetic molecules a representation of 
the evolution by a Markov process achieves the required probability model that is used to study the 
probability density function (PDF) of the order parameter, i.e. the magnetization. The existence 
of one or more modes in this PDF is an indication for the superparamagnetic transition of the 
cluster. This allows us to determine the factors that influence the blocking temperature, i.e. the 
temperature related to the change of the number of modes in the density. It turns out that for our 
model, rather than the evolution of the system implied by the Hamiltonian, the high temperature 
density of the magnetization is the important factor for the temperature of the transition. We find 
that an initial probability density function with a high entropy leads to a magnetic cluster with a 
high blocking temperature. 
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I. INTRODUCTION 



Heating a ferromagnet will make it paramagnetic at a sufficiently high temperature. By 
cooling and magnetizing, it can be made ferromagnetic again. This simple fact of everyday 
magnetism is also true on a nano-scale. As a matter of fact the existence of this phenomenon 
has been predicted by Neel [l| and is indicated by superparamagnetism. The temperature 
of the transition is the blocking' temperature and is experimentally known in molecular 
magnets and magnetic clusters . It casts a shadow on possible applications 

because this temperature is very low for most materials. A theoretical estimation of this 
temperature could offer a better understanding of the factors that have an influence on its 
value and would indicate a way to devise clusters with an appropriate blocking temperature. 

As far as we know, there is no well established procedure to obtain this information, 
especially when the system is in the quantum regime, statistically non-extensive and the 
blocking temperature has to be determined from some kind of a phase-transition in a finite 
system. Most theoretical approaches are based on the Fokker-Planck equation or make use 
of Langevin dynamics t 8j. These methods assume that the magnetization is a continuous 
variable, an assumption that is inappropriate for systems whose magnetic degrees of freedom 
are small. Furthermore for molecular magnets embedded in a crystalline structure the 
iamiltonian H describing the spectral properties of the basic molecular unit is well known 
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12j. It is used to relate the spectral properties of the system to a lot of observed 

fl nnn 



responses, e.g. quantum tunneling of magnetization |4|, |ll|, |l3j, ll4( . Therefore a study of the 
blocking temperature should take the information obtained from the spectrum into account. 

Recently some substantial progress was made in non-extensive statistical physics for nu- 
clear clusters In the same line of thought and more important for the approach we 
will use here, is the observation by Gulminelli 16[ about the equivalence of the location of 
the zeroes of the partition function in the complex temperature plane and the modes of the 
probability density function (PDF) of the order parameter. Investigations along this line 



include: metallic clusters |Tj|, pairing in nuclei 



capacity jsoL It was first shown that for non-extensive systems the location of Yang 



^ HQ, 



and clusters with a negative heat 
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24j . Subsequently the 



Lee zeroes reveals the transition from one phase to another 
equivalence between the transition and the modality of the PDF of the order parameter 
offers a very direct way to obtain the temperature range where one may expect a change 
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from ferromagnetic to paramagnetic behavior in the case of the nano-magnets. Indeed ferro- 
magnetic behavior corresponds with a two-mode PDF for the magnetization, paramagnetic 
behavior with a one-mode PDF. 

In order to calculate the PDF of the magnetization for a molecular magnet or cluster 
we need a probability model able to describe cooling and heating and compatible with that 
Hamiltonian. Our starting point is a generalized propagator of the form: 

K(m,Tn',/3) = (m \ exp(—j3H) | m')—^. (1) 

9m' 

In some cases I25I l^. the appropriate choice of the states m and the function g m will lead 
to an evolution equation for the transition probability of a Markov chain. For the given H 
we could indeed map the Hamiltonian description to a Markov process in a space spanned 
by the possible outcomes of the magnetization of the cluster. The process, only defined for 
(3 > 0, describes the cooling of the system. In order to obtain the inverse evolution, the 
Bayes theorem can be invoked to define the adjoined heating process. In our example the 
dimensionality of the process remains low. Therefore only standard numerical methods are 
necessary to perform the calculations. 

Our task will be to obtain this PDF for a system with an evolution described by the 
process generated by the appropriate Hamiltonian. Once this calculation is done, it is 
our hope to find a functional relation between the parameters in the Hamiltonian and the 
blocking temperature. It turns out that other factors play a more important role. In 
order to identify these factors we will introduce in the second section, the Hamiltonian, the 
mapping to a Markov process, how to describe cooling and heating and outline the numerical 
techniques to obtain the PDF of the magnetization. In the third section we present a study 
of the PDF of the magnetization. We point out, by means of several cases, that preparation, 
exactly as with ordinary magnets, is of fundamental importance. We summarize our results. 
In the last section we discuss the method, the results and try to give some perspectives. 



II. CONSTRUCTION OF A PROBABILITY MODEL 



The model Hamiltonian for a molecular magnet takes the standard form 
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H = D(S Z ) 2 + E((S X ) 2 - (S y ) 2 ) + g^ B (b z S z + b x S x ). (2) 
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where D and E are the diagonal and non-diagonal anisotropy constants, respectively and 
b x (b z ) is the external field applied along the specified axis. Rewriting (J2J) and replacing the 
cartesian spin components by raising and lowering operators, one can split the Hamiltonian 
into three parts: 

H = H + V 1 + V 2 (3) 

The term Hq is diagonal in the z-component of the spin. The term V\ of the Hamiltonian is 
proportional to the x-component of the magnetic field and linear in the lowering and raising 
operators: 

Vi = gj^ (s + +s-) . (4) 

The term V 2 is proportional to E and quadratic in the lowering and raising operators: 

V 2 = j (S+S++S-S-) . (5) 

For some molecular magnets this term is forbidden by symmetry and a fourth order term is 
introduced. Our model allows that a term containing a product of a lowering and a raising 
operator can be rewritten in terms of the z-component. 



A. The evolution equation 

In order to derive the evolution equation let us consider first V\. The states are determined 
by two numbers N and m: in terms of a spin | system N represents the number of spins, the 
number m is an integer of the interval [0, N] and represents the number of spins with their 
magnetic component in the opposite direction. In general the magnetic quantum number is 
given by: 

( N 

S z \m) — m 

We have chosen the number m to represent the states because it allows a direct way to mon- 
itor the transitions by the counting the number of spins in the opposite direction. Lowering 
or raising the number m is achieved by the other spin operators that act as follows on a 
state: 

S + \m) = y/(N -m)(m + 1) \m + 1) , (7) 
S~ \m) = a/ (N — m + 1) m \m — 1) . (8) 
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The matrix elements of the evolution equation in terms of these states are readily obtained 
from: 

(k\ V x \m) = (y(N-m) (m + l)4 m+1 + y/(N-m + l)l5 k ro _x) . 

The square roots in this representation arise through the normalization of the states neces- 
sary for the Hermitian character of the Hamiltonian and they are eliminated by the choice 
g m = a/ (2s — m) \ (m)! in (JTJ meaning that the normalization factor is transferred to the ket 
of the quantum state. 

The matrix elements of Vi are straightforward to calculate using the same technique and 
with the same g m the propagator becomes 



/ x , / T n , a f m\(N -m)\ I 2 

K(m,m',r) = (m | exp(-riJ) | m') | ^ _ J } , j . (9) 

and the evolution equation of the model reading now: 

dpK(m,m! ', (3) = w(m,m + 8)K(m + 5,m' , (3) — e m K(m,m' , j3), (10) 

<5=±1±2 

does not contain square roots anymore. The energy spectrum of the diagonal part of the 
Hamiltonian (0) is given by 

n .iV-2m. 2 N-2m 
e m = D{ ) 2 -gfi B b z (11) 

and the rates w(m, m + 5) are derived from the term V1 + V2. These rates are: 

w(m, m + 1) = b(N — m) (12) 

w(m, m — 1) = bm 

w(m,m + 2) = E(N — m)(N — m — 1) 

w(m,m — 2) = Em(m — 1) 

where b = gfisbx and E are chosen to be positive. For our model (j2J) it can be done without 
loss of generality. It should be noted that the transitions do not lift the degeneracy of a 
configuration with magnetization characterized by m- |. The evolution equation (jl(J|) can 
be written as a matrix equation with dimension (N + 1) x (N + 1). The formal solution of 
this equation is: 

K(/?)=exp(-/?(E-W)), (13) 



where E is a diagonal matrix with e m on the diagonal. The matrix W has the rates w(m, I) 
as elements. 

Although for modest values of N the matrix (J 13)) can be calculated directly using symbolic 
algebra or numerically using standard techniques, the calculation has to fit in the probability 
theoretical setting mentioned in the introduction. This can be done by calculating the 
energy and multiplicity of the quantum states and use the Gibbs rule to relate these energy 



eigenvalues to probabilities. This road has been followed in the references 
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23, 



30. 31. 



The alternative is to introduce a Markov representation for the quantum evolution 



(in imaginary time or equivalently in inverse temperature) of the system. 



B. The Markov representation 

Before we illustrate the method used to transform the system described by |T0|) . we will 
argue on the opportunity to choose a Markov chain for the evolution. First of all the form 
of the Hamiltonian allows this kind of representation, as far as we know the mapping is 
limited by constraints on the transition rates and therefore not general. Furthermore once 
a mapping can be found it is straightforward to construct a probability model without 
further assumptions. In other words the mapping makes the choice of eigenstates of the 
magnetization as the state space optimal with respect to the evolution inducing transitions 
in the state space. The Markov property assures thereby that transition probability to go 
to the next state is independent of the history preceding the transition into the actual state. 

The evolution equation can be split into two parts: 

d T K(m,m' , r) = w(m,m + 5)K(m + 5,m\r) — u m K(m,m' ,r) — (e m — u m )K(m,m' ,r), 

<5=±1±2 

(14) 

where the following part describes a Markov chain 

d T K°(m,m' ,r) = w(m, m + S)K°(m + 5, m', r) — v m K°(m, m', r) (15) 

5=±1±2 

provided v m is given by: 

v m = w(m,m + 6). (16) 

<5=±1±2 

On the condition that (e m — u m ) > it is possible to use K(m, m', r) to define a probability 
model for the evolution of the system. Under this condition it is straightforward to show 
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that K(m,m' , r) > 0. A heuristic proof is given in the appendix. However in order to 
fall in the category of Markov chains it is necessary that Estates K(m, state, r) = 1. For 
K°(m,m' , r) one can easily prove this property. 

For K(m, m', r) one can construct an equivalent Markov chain by adding an extra state 
c to the state space. Using this state which is an absorption state sometimes called the 
coffin-state one obtains: 

^K{m,m',T) + K{m,c,r) = 1 (17) 

m' 

and defining K(c, m', r) = for all ml ^ c and K(c, c, r) = 1 one obtains the new process. 
This means that we obtained a probability model with the conditional probability given 

by: 

P[M T = m | M T i = to'] = K(m', m,r - r'), (18) 

provided r > r' . M T is a random variable that takes [to = 0..N] or c as a value. The 
conditional probability is derived from the solution of the evolution equation. In terms 
of the inverse temperature it should be noted that the evolution is only represented by a 
Markov process when the increments are positive. 



C. Cooling and heating 

Cooling is described in a straightforward way: noting that in a model that describes the 
evolution to equilibrium the variable r represents the inverse temperature, a solution of the 
evolution equation allows to calculate the conditional probability to reach a given state at 
t given the state of the system at r' . Because r > r' one deals with cooling. 

In order to describe heating however a further elaboration of the probability model is 
necessary. Assume that the inverse temperature f3 is larger than /3' and suppose we know 
that the system is in a given state at j3 then we may ask ourselves: what is the probability 
to find it in a state to' at (3'1 The quantity to obtain is: P[Mpi = m! \ Mp = m] . Using 
Bayes' theorem one finds: 

P[M P , =m'\Mp = m}= P[M P — to | Mpr — ^' fp^ Z^ ( 19 ) 
Noting that if we know the density of the states at r = the probabilities P\Mp = to'] and 
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P{M[3 = to] can be calculated. They are given by: 

P[Mf3 = to] = p i M l3 = to | M = k]P[M = k], (20) 

k 

and the transition probability for heating can also be obtained from the solution of evolution 
equation provided the initial PDF is known. 

In this section we showed that under quite general assumptions about the parameters of 
the model the evolution of the model under cooling or heating can be cast into a Markov 
chain. 

III. THE PDF OF THE MAGNETIZATION 

A key ingredient in our probability model is the density at (3 — 0. In order to obtain the 
PDF of the magnetization at a given temperature either by heating or cooling the system, 
one has to know the initial PDF for the magnetization: P[M = k). Before we analyze the 
temperature dependence of these initial PDF's we will comment on two possible choices. 
Both choices are based on the principle of equal "a priori" probabilities. If the system can 
be described by a single spin statistically, the probability density is uniform for all states 
that are participating in the evolution: 

P[Mo = k] = ] Jp I . (21) 

If the system allows different spin configurations leading to same magnetization, one may 
assume that the M = is realized by all spins in the same direction and the M = N by all 
spins in the opposite direction. The other states are combinations of (N — k) states with spin 
up and k states with spin down. The initial density is calculated by considering S? = ±| 
spins with S z = J2i The total number of spins pointing up determines the value of the 
total S z uniquely. The number of ways to realize this number is easily obtained. Taking to 
into account that for a spin | the probability that a single spin Sf is pointing up, is |, the 
binomial density follows: 

P[M = k] = 2- N (^j. (22) 

These choices correspond respectively with a single large spin or with many spins all having 
an evolution described by the same Hamiltonian. In order to distinguish between the choice 
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of the initial density an index will be introduced: we will use LS for an initial uniform 
density (|21|) and MS for the binomial density (|22p. 

The calculation of the PDF of the magnetization is straightforward: first the matrix ele- 
ments of (|13|) are calculated by a direct method. Once these matrix elements are known the 
transition probabilities of the Markov representation (|18|) are obtained, taking the relation 
between of S z and Mp determined by their definition © into account. 

Here we will describe the cooling and heating of the LS model first remarking that for 
b z = and realistic values of D, E and b x there is no blocking temperature i.e the PDF 
remains bimodal ((TJ. We note further that from the first and second moment of the PDF, 
i.e. the mean magnetization and specific heat, it is almost impossible to infer that these are 
moments of a bimodal density. 
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FIG. 1: The probability density of the magnetization for the single spin model The 

single spin model has N + 1 realizations of the magnetization, each of these realizations is equally 
probable at high temperature. Using this as initial condition (|21j) the density of the magnetization 
at a given temperature can be obtained. As illustrated here for three typical temperatures, it is 
seen that the density remains bimodal in the whole temperature range. 

Next we will describe the heating and cooling of the MS model and see that this model 
has a blocking temperature: changing the bimodal structure of the PDF into a PDF with 
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FIG. 2: The probability density of the magnetization for the many spin model The many 
spin model has 2 N realizations of the system leading to a binomial density of the magnetization at 
high temperature. Using this as initial condition (|22|) the density of the magnetization at a given 
temperature can be obtained. As illustrated here for three typical temperatures, it is seen that the 
density is bimodal in the low temperature range but has a single mode in the high temperature 
range. The change from a density with a single mode to a bimodal density occurs at the blocking 
temperature. 

one mode on heating from the low temperature phase and changing the PDF with one mode 
to a bimodal PDF on cooling from the high temperature phase (figJ2J). Two remarks are 
in order now: firstly we note that the blocking temperature has the same numerical value 
for cooling and for heating within the accuracy of our calculation, secondly the order of the 
transition based on the continuity of a distribution function cannot be used to classify the 
transition because the model has a discrete distribution function. 

From the preceding analysis we arrive at a provisional conclusion: the blocking tempera- 
ture of the model crucially depends on the initial density. This density is obtained by using 
the principle of equal a priori probabilities for all the states participating in the evolution 
of the system. In the sense that we have to know the multiplicity of a state with M = m 
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in order to assign an initial PDF. 
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FIG. 3: The initial probability density of the magnetization As the spin number £ decreases 
the change of the initial magnetization density from a binomial distribution £ = 20 to an uniform 
distribution £ = (|23|) is illustrated for N = 20 spin cluster. The initial entropy values and £ 
characterizing the initial distributions are mentioned in the legend. 

In most molecular magnets or clusters different building blocks can be distinguished: 
this has led us to consider a third model where we assumed that the magnetization can 
be decomposed into a sum of two independent components. The first part initially has a 
binomial density: 

M = M £ + Mg~ l 
where i is the number of spins. The initial density of the first part is: 

P[M* = *] = 2-<Q, 
the density of the second part is taken from a single spin and is uniform 

p[m^- £ = i] = i/(N-e + i). 

This leads to an initial density for the magnetization given by a convolution: 

P[M = m] = e(m = I + fc) (jV 2 _^+i) - ( 23 ) 
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Once the initial density is chosen (figEJ) it is straightforward to calculate the PDF of 
the magnetization and to determine the blocking temperature. It should be noted that the 
new model interpolates between the two models previously introduced. Taking i = is 
the large spin model and taking I = N gives the many spin model. In order to relate the 
blocking temperature of the model and a characteristic of the initial distribution, we have 
calculated the initial entropy and we have plotted the obtained blocking temperatures (figEJ) 
in function of the initial entropies. 



0.4 




FIG. 4: The blocking temperature versus the initial entropy The blocking temperature 
is derived from the change of the number of modes in the density of the order parameter. This 
change occurs at a different value of (3 for a different initial density (see figEJ)- The initial density 
is a convolution of a uniform density with a binomial one, characterized by S, the entropy of the 
initial density. 



The same information can be given in a bifurcation diagram (figEJ) following the modes 
of the density. The blocking temperature of the preceding plot (figEJ) is obtained by deter- 
mining the value of (3 where the bifurcation occurs. 



IV. DISCUSSION AND CONCLUSION 



In section 2 we mapped the evolution described by the Hamiltonian of the magnetic 
cluster on a continuous time Markov chain, time refers here to the ordering parameter /3, 
the inverse temperature. This mapping is a very appealing way to obtain a probability model 
compatible with the Hamiltonian of the system. This model is a necessary requirement in 
order to use the analysis of phase transitions or finite systems mentioned in the introduction. 
In order to have a complete description, not only cooling but also heating has to be considered 
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FIG. 5: The bifurcation of the modes for the density of the order parameter The modes 
of the density of the magnetization are plotted in function of the inverse temperature j3 for different 
initial densities (J23J) characterized by £. 



as a process. This is achieved by introducing an absorbing state, the so-called coffin state. 
Once this formalism is chosen, the procedure to calculate a PDF at a given (3 requires the 
knowledge of an initial PDF. The reason for this requirement comes from the formulation of 
the evolution in terms of transition probabilities, conditional by definition. We showed that 
for different choices of the initial PDF but with the same set of transition probabilities the 
blocking temperature critically depends on the initial condition. Our conclusion is reached 
by the use of a mathematically well-established procedure for Markov systems. We have to 
mention a possible alternative and eventually discuss the differences if any in the conclusion. 

The more traditional approach would be to calculate the trace of the propagator to obtain 
the partition function and to use the diagonal elements of the propagator explicitly in the 
Gibbs assumption to obtain a PDF: 

^ r ,^ , K(m,m,3) . . 

P[ Mf3 = m]= u m ■ (24) 



We have to introduce here v m |35[, the multiplicity of the state m. In the present formalism, 
the initial PDF contains this information. Comparing ()24|) with the PDF calculated in the 
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FIG. 6: Comparison of PDFs for Gibbs assumption and our probability model The 

numerical difference between the PDF obtained in the Gibbs assumption (□) (|24jl and in our 
probability model(T) (|2"U|) introduced in section 2 at a given inverse temperature (3. 



preceding section, we find that they look very similar but there are numerical differences 
as can be seen in fig©. This means that also the blocking temperature changes slightly 
passing from one approach to the other. In view of the crucial factor played by the initial 
condition for the PDF and its relation to the multiplicity of the states, it should be noted 
that the initial condition is inherently related to the use of a Green function, i.e. the 
propagator, while in the Gibbs assumption (|24|) it is implicit in the summation over the 
states. A profound analysis of the origin of the difference between (J2IJ) and (f2T?j) would 
require additional research. The main difference in approach is the choice of the states: if 
the eigenstates of the Hamiltonian are considered as states of the probability model there is 
no reason to construct a Markov representation and the invariance of the trace with respect 
to cyclic permutations can be used to avoid explicitly a diagonalisation. The probability 
model that we considered is formulated in terms of the eigenstates of the order parameter, 
i.e. the eigenstates of S z . Now S z and the Hamiltonian do not commute, therefore there is no 
common base that can serve as states of a system where both operators are diagonal. From 
this point of view generating the evolution of the system in terms of the states characterised 
by m cannot avoid to take the transitions into account. 

The model for a cluster containing for example N spin | is clearly a many spin model, 
however most molecular magnets are described by a theoretical model that relies on a single 
large spin description. Fortunately both models obey the same dynamics. They have the 
same spectrum but a different dimensionality leading to a different blocking temperature. 
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In the many spin model with dimension i the one-to-one correspondence between the mag- 
netization and the states implies that a specific state characterized by m spins with a s z 
component down, have a multiplicity given by ( ) in order to conserve the total number of 
states 2 . As indicated in our investigation about phase transitions in systems with finite 
degrees of freedom, the multiplicity of the states is an essential ingredient. Compare for 
instance the magnetization of the single large spin model where the spin states have mul- 
tiplicity with the many spin model with a binomial density as initial PDF. It is clear that 
the binomial multiplicity can be related to dimension. Take a harmonic oscillator in one 
dimension, a state with the energy level nftw has the multiplicity 1, while the multiplicity 
of that state is a function of n in higher dimensions. This suggests that the single large 
spin model can be seen as a one-dimensional realization of the system. According to the 
Van Hove theorem, it should not have a phase transition and indeed there is no transition. 
However it was a surprise, for us, that the equilibrium probability density of the magnetiza- 
tion resembles that of the low temperature phase of the higher dimensional model, strongly 
indicating that one-dimensional models are low-temperature by nature. 

To summarize we do not change the evolution of the system, we change only the initial 
density of the magnetization by making assumptions on the magnetic structure of the cluster. 
The simplest magnetic structure in terms of degrees of freedom - each state can be realized 
only in one way - leads immediately to a bimodal density and therefore there is no transition 
at finite temperature. The structure with a large degree of freedom - the total number of 
realizations is 2 N - has a binomial initial density and leads to a finite blocking temperature. 
Interpolating structures with a density that arises from a convolution of binomial densities 
and uniform leads to larger blocking temperatures. This means that the model indicating 
the smallest blocking temperature is attributed to the model with the largest number of 
realizations of initial states. Less realizations of the initial states leads to larger blocking 
temperature. If this trend remains true for more realistic descriptions of magnetic clusters, 
it can be an important factor to take into account in the design of the cluster. 

For a finite system with many degrees of freedom but a relatively simple evolution we 
could calculate the blocking temperature by mapping the system on a Markov chain in 
order to obtain a probability model with a change in the number of modes of the probability 
density function of the order parameter. This approach rests on the construction of such 
a probability model. It is clear that new methods are required in order to achieve such a 
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mapping for systems described by a more realistic Hamiltonian. It should be noted that 
simulating the cluster belongs to the same category. In the sense that: it is relatively easy 
;o take the diagonal elements of the propagator into account using a Metropolis algorithm 
3f|, but combining this algorithm with a Markov chain incorporating transitions in the 
simulation requires further investigation. 

APPENDIX A 

In order to show that the diagonal elements of K can be used to define a probability, the 
trace should exist and the elements should be positive. The following evolution equation is 
considered first: 

d T K (m, m', r) = w(m, m + 5)K°(m + 5, m', r) — w(m,m + 5)K°(m,m',r). 

<5=±1±2 <5=±1±2 

(Al) 

This is the evolution equation of a Markov process defined on the states m £ [0, N] provided 
that the w(m,m + 5) are positive. This means that with the appropriate normalization 
K°(m, m', r) is a transition probability hence positive and finite. Using the transition prob- 

that the solution of (fTT)|) is given by an expectation 



ability K°(m, m', r) it can be shown 



of a Feynman-Kac functional for V(m) = e m — Yls=±i±2 w i m i tu + S) 

K(m,m',T)=S m [I{M{r)-m')exp{- f V(M(s))ds)} (A2) 

Jo 

where M(t) is the random process of a Markov chain with K°(m,m' , r) as a transition 
probability and £ m is the expectation over that process. The random variable I(M(t) — m!) 
is an indicator necessary to satisfy the condition that the system is initially in the state m' . 
This shows that K(m,mf , r) can be decomposed into a sum of positive terms. 
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